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The results of muon-spin relaxation and heat capacity measurements on two pyroxene compounds 
LiFeSi2C>6 and NaFeSi2 0e demonstrate that despite their underlying structural similarity the mag- 
netic ordering is considerably different. In LiFeSi2 06 a single muon precession frequency is observed 
below Tn, consistent with a single peak at Tn in the heat capacity and a commensurate magnetic 
structure. In applied magnetic fields the heat capacity peak splits in two. In contrast, for natural 
NaFeSi2 0e, where multiferroicity has been observed in zero- magnetic-field, a rapid Gaussian depo- 
larization is observed showing that the magnetic structure is more complex. Synthetic NaFeSi2 06 
shows a single muon precession frequency but with a far larger damping rate than in the lithium 
compound. Heat capacity measurements reproduce the phase diagrams previously derived from 
other techniques and demonstrate that the magnetic entropy is mostly associated with the build up 
of correlations in the quasi-one-dimensional Fe 3+ chains. 



PACS numbers: 76.75.-H, 75.50.Ee, 75.85.+t 

I. INTRODUCTION 

Multiferroic materials demonstrating coupled mag- 
netic and ferroelectric order have once again become an 
active field of research, since they offer both interesting 
physical properties and the possibility of technological 
applications. 1-6 While an increasing number of multifer- 
roic materials have been discovered in recent years^— 
and much progress has been made in finding general rules 
to describe the origins of this effec t 14 ' 15 it is not always 
possible to predict if a given material will be multifer- 
roic. Isostructural series have already provided consider- 
able insights into multiferroicity, notable examples being 
the hexagonal and orthorhombic manganites i 16 ' 17 Com- 
peting magnetic interactions and a strong magnetoelas- 
tic coupling are both known to favour multiferroicity. In 
this context the discoveries of multiferroicity in the py- 
roxene compound NaFeSi2C>6 and magnetoelectricity in 
LiFeSi2C>6 have suggested that this geologically common 
family may offer more multiferroic compounds, as well as 
providing an opportunity to study isostructural materials 
with different spins and magnetic exchange constants^ 

Pyroxene compounds have chemical formulae 
A+M 3+ (Si, Ge) 2 6 and chains of M 3+ ions surrounded 
by oxygen octahedra lie along the crystallographic 
c-axisi 1 ^ The M 3+ chains are connected by (Si, Ge)C>4 
tetrahedra. This structure is shown in Fig. [T] Most 
magnetic members of this family show Neel ordering at 
low temperature, a notable exception being the orbitally 
assisted spin-Peierls transition seen in NaTiSi206' 2 ° l21 
Common to both the Neel ordered and spin-gapped com- 
pounds is the dominant intrachain exchange interaction 
giving quasi-one-dimensional magnetic properties. 




FIG. 1: (Color online) Structure of LiFeSi206 showing the 
Fe 3+ chains running along the c-axis linked by Si04 tetrahe- 
dra. The structural data come from Ref. |22j. 



LiFeSi206 has a Neel temperature of 18 K and there 
is no pyroelectric current without an applied magnetic 
field. 18 Applying a magnetic field along the c-axis re- 
duces the temperature of the peak of the magnetic sus- 
ceptibility to 14 K at 14 T and measurements of the 
pyroelectric current, lb, show a peak which follows the 
same magnetic field dependence as that in the magnetic 
susceptibility. Smaller peaks in lb at higher tempera- 
ture were also observed but their origin is unclear.^ The 
magnetic structure has been determined to be antiferro- 
magnetically coupled ferromagnetic chains with magnetic 
space group P2i/c'M This magnetic structure allows for 
magnetoelectric effects consistent with those observed. 
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The isostructural compound LiCrSi2C>6 was found to 
have comparable magnetic and magnetoelectric proper- 



ties 



18.23 



The situation in NaFeSi206 is rather more complex, 
largely because of the differences observed between nat- 
ural and synthetic samples. Natural samples, which are 
known to contain impurities, show two phase transitions 
in zero magnetic field: at 8 K to a collinear magnetic 
structure and at 6 K to a ferroelectric (P \\ b) phase. In 
fields above 4 T a ferroelectric (P || c) phase was ob- 
served below 5 KJ£ Synthetic samples showed a similar 
magnetic structure to LiFeSi2 0g, with antiferromagneti- 
cally coupled ferromagnetic chains, albeit with evidence 
for a further incommensurate modulation to this struc- 
ture that could not be determined^ Given that it was 
not possible to index all the magnetic Bragg peaks, the 
results of ab initio calculations, and the multifcrroicity 
observed in the natural samples it seems likely that the 
magnetic ordering is actually helical and incommensu- 
rate^ 

More detailed ab initio calculations for a broad range 
of pyroxene compounds, including those we study here, 
were carried out by Streltsov and Khomskii— . They 
modelled the exchange constants in terms of an intra- 
chain exchange J, and two interchain exchange constants, 
J\ and J2 , all of which were found to be antiferromagnetic 
for both compounds. The calculations suggest J Ll = 7 K, 
J x Li = 1.9 K, and J 2 Li = 3.4 K; and J Na = 8.5 K, 
Jf a = 0.8 K, and J 2 Na = 1.6 K.^ These values suggest 
that the magnetism in LiFeSi206 is likely to be more 
three-dimensional than that in NaFeSi206 and the dif- 
ferent exchange constants may have an even more signif- 
icant effect on the fine details of the magnetic structure 
and any magnetostriction. 

In this paper we investigate synthetic samples of the 
two pyroxene compounds LiFeSi2 06 and NaFeSi2 06, and 
a natural sample of NaFcSi206, using heat capacity and 
muon-spin relaxation measurements. These probe the 
change in magnetic entropy around the phase transitions 
and the local magnetic field distributions within the sam- 
ples. While there are some underlying similarities in the 
magnetic properties, the effects of the changing exchange 
constants and the presence of impurity-induced disor- 
der in the natural sample are clearly evident in the data 
recorded by both techniques. 
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FIG. 2: (Color online) Heat capacity measurements on 
(a) LiFeSi2 06, (b) synthetic NaFeSi2 0e, and (c) natural 
NaFeSi2 0s. The lines in the main panels show the fitted 
lattice terms described in the text and the insets show the 
variation with applied magnetic field close to the magnetic 
ordering transitions described in the text. 



was obtained by crystallisation of glassy NaFeSi2C>6 that 
was prepared using high-temperature flux. 



II. EXPERIMENTAL 

A. Samples 

Our natural sample of NaFeSi2 06 was cut from 
the same crystal that was used in Ref. [l|| Elec- 
tron microprobe analysis has shown that the composi- 
tion is Nai.o4Feo.83Cao.o4Mn .o2Alo.oiTio.o8Si20 6 i^ The 
synthetic sample of LiFcSi2C>6 was composed of small 
translucent single crystals grown from melt solution, see 
also Ref-ITsl The powder sample of synthetic NaFeSi2C>6 



B. Heat capacity measurements 

Heat capacity measurements were made using a Quan- 
tum Design Physical Properties Measurement System 
(PPMS), employing the two-tau relaxation method, in 
magnetic fields between and 10 T. The data for both 
compounds are shown in Fig. [2J They show clear peaks 
associated with the magnetic ordering transitions found 
using other techniques^ In compounds such as these, a 
(partial) hump in the heat capacity is observed above 
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Tn as the correlations build up, and then a peak is 
superimposed upon this as the system enters a three- 
dimensionally ordered state. Some information concern- 
ing the dimensionality of the system can be obtained from 
the form of the hump and the relative size of the peak3& 
The lattice contribution to the heat capacity of each com- 
pound was modelled using one Debye and two Einstein 
components. Parameters derived from fitting this form 
to the data above 50 K are given in Table HI 

In zero applied field LiFeSi20g shows the build up of 
short-ranged correlations in the chains from well above 
Tn and a single peak in the heat capacity at Tn = 18 K. 
Having subtracted the lattice contribution, we estimate 
that the integrated magnetic entropy up to 50 K is 
11.2 Jmol K , with around 80 % accounted for by 
the short-range correlations. The heat capacity mea- 
surements on LiFeSi2 06 show an unusual hysteresis with 
a magnetic field applied perpendicular to the ab plane. 
Measuring in successively increasing magnetic fields up to 
10 T repeated the peak in the data seen in zero-field [0 T 
(A)] within the experimental error. However, measure- 
ments in successively decreasing fields, while cooling the 
sample from around 2 x Tn, gave the two field dependent 
steps in the heat capacity shown in the inset to Fig. [5Ja) . 
No pattern is evident in different positions of these steps 
in different fields. The two features were accompanied by 
small amounts of latent heat (evident in the poorer fits to 
the raw thermal relaxation data recorded by the PPMS) 
and persist down to zero applied field [0 T (B)]. This 
behavior suggests short range order persists well above 
Tn and produces hysteresis in the sample when fields are 
applied. 

The heat capacity data on synthetic NaFeSi20g are 
shown in Fig. [2fb) and take a similar form to the data 
for LiFeSi2 06, with a significant magnetic heat capacity 
well above Tn- A clear peak is found at 6.6 K, rather 
broader than in LiFeSi2 0g and a little lower than the 
transition temperature found in the /iSR measurements. 
The peak broadens with increasing field but does not 
move significantly. Our data are in excellent quantitative 
agreement with those reported previously by Ko et al.— 
for their measurements on a synthetic sample. The inte- 
grated magnetic entropy up to 50 K is 11.2 Jmol K -1 
with around two-thirds of this appearing to be associated 
with the build up of correlations within the chains. 

In natural NaFeSi2C>6 [Fig. Efc)] the field- dependent 
part of the heat capacity forms a much smaller fraction 
of the feature around 7 K, which is dominated by the 
build up of correlations in the chains. Examining the 
data closely shows two small peaks at 8 K and 6 K corre- 
sponding to the magnetic transitions and these merge as 
the magnetic field is increased, consistent with the phase 
diagram proposed in Ref. [Lsl . The magnetic entropy in- 
tegrated up to 50 K is approximately 13.9 Jmol _1 K _1 , 
approximately 90 % of which is accounted for by short- 
ranged ordering. It is notable that the field dependent 
heat capacity of LiFeSi2 06 shows two peak features in 
applied field that are similar to those seen in natural 



Sample 


LiFeSi 2 6 
synthetic 


NaFeSi 2 6 
synthetic 


NaFeSi 2 6 
natural 


0D (K) 
fcsi (K) 

0E2 (K) 


621(18) 
192(3) 
1210(50) 


370(20) 
190(10) 
700(50) 


530(20) 
165(5) 
1050(50) 



TABLE I: Parameters and their statistical errors derived from 
fitting one Debye and two Einstein components to the heat 
capacity data above 50 K. 



NaFeSi2C>6. The origin of the two features is at present 
unknown and neutron scattering experiments in applied 
field should be carried out on LiFeSi2C>6 to discover 
whether a magnetic transition is involved. 



C. fiSR. measurements 

Our positive muon-spin relaxation (/iSR) measure- 
ments^ (r M = 2.2 (is, 7 M = 2tt x 135.5 MHzT -1 ) were 
carried out on the General Purpose Surface-Muon In- 
strument (GPS) at the Paul Scherrer Institute, Switzer- 
land. Samples were mounted on a low background sam- 
ple holder with aluminized mylar tape to minimize the 
background from muons stopping outside the sample. To 
measure the time evolution of the muon spin polariza- 
tion, emitted decay positrons were collected in detec- 
tors placed forward (F) and backward (B) relative to 
the initial muon spin direction (antiparallel to the beam 
momentum). The muon decay asymmetry is defined in 
terms of the count rates in the two detectors (Ap and 
A B ) as: 



A(t) 



N-p(t) - aN B (t) 
N F (t)+aN B (t)' 



(1) 



where a is an experimental calibration constant related 
to the relative efficiency of the detectors. 

The muon spins are sensitive to both static and fluc- 
tuating local magnetic fields at their stopping positions 
inside the material, and these affect how the form of the 
muon decay asymmetry changes with time. In the para- 
magnetic phase of each compound the muon relaxation 
is well described by a single exponential relaxation. In 
LiFeSi 2 6 and synthetic NaFeSi 2 6 [Fig. [3] (a) and (b)] 
we observe coherent muon precession below Tn consistent 
with long range magnetic order and quasistatic magnetic 
fields at the muon stopping site. The data are well de- 
scribed by the function: 



A(t) = A x e~ Xlt cos(27rz/t) + A 2 e 



-\ 2 t 



(2) 



The first term describes damped muon precession around 
quasistatic local fields (B = 27w/7 /J ) perpendicular to 
the muon spin polarization and the second term is an 
exponential relaxation, of rate A2, due to fluctuations 
flipping the spins of muons having a non-zero spin com- 
ponent along the local magnetic field direction. The val- 
ues of v derived from fitting the asymmetry data are 
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FIG. 3: (Color online) Muon asymmetry data above and be- 
low the magnetic transitions for: (a) LiFeSi2C>6, (b) synthetic 
NaFeSi206, and (c) natural NaFeSi2C>6. Below Tn we fit with 
Eq. [2] for (a) and (b), and Eq. g] for (c). Above T N the re- 
laxation is exponential. The width of the time bins in the 
asymmetry histograms have been increased for clarity. For 
the synthetic NaFeSi206 measurements were done with the 
initial muon spin rotated differently relative to the detectors, 
leading to the lower asymmetry values. 



shown in Fig. |U[a). We find that about 2/3 of the asym- 
metry is associated with the oscillating signal, consis- 
tent with the fact that in a polycrystalline sample 2 /3 of 
the muons will find local magnetic fields perpendicular 
to their spin polarization and 1/3 will experience fields 
along their spin direction that can only lead to depolar- 
ization if fluctuations are present. In LiFeSi2 0g Ai and 
A2 are almost temperature-independent. The parameters 
extracted from the /iSR data analysis are presented in 
Fig. [4] As shown in Fig.|4ja), the precession frequencies 
in LiFeSi2 06 and synthetic NaFeSi2C>6 are well described 



by the phenomenological function: 

v{T)=v{Q){l-{T/T^) a f. (3) 

For LiFeSi 2 6 , i/(Q) = 28.9(3) MHz, T N = 18.50(1) K, 
a = 1.6(2), and (3 = 0.26(2). This implies that the T ->• 
internal field at the muon site is approximately 0.2 T. For 
synthetic NaFeSi2C>6 the frequency is less well defined 
because the oscillations are far more strongly damped 
and, constraining a to the value found for the Li sample, 
we can fit v{0) = 27(1) MHz, T N = 7.07(5) K, and /3 = 
0.28(4). 

Data for natural NaFeSi2C>6 takes a different form to 
that for the other two samples [see Fig. [3] (c)]. Here 
we see no coherent muon precession, suggesting a large 
range of quasistatic fields occurs at muon stopping sites. 
A random distribution of static magnetic fields leads to 
a Kubo-Toyabe function^, which shows a dip and re- 
covery in the asymmetry. Even with the addition of a 
damping term it is not possible to get reliable fits to the 
measured data using such a fitting function. Instead, we 
can effectively describe the measured asymmetry using 
a rapid Gaussian relaxation to describe the effect of the 
quasistatic fields and a slow exponential that describes 
the 1/3 tail expected for the Kubo-Toyabe function: 

A{t) = A x e-^ e +A 2 e- Xt . (4) 

In analogy with Eq.[2]the first term describes the incoher- 
ent precession about large static magnetic fields and the 
second term describes spin-flipping of muons with their 
spin direction aligned along the local magnetic field. The 
experiment on natural NaFeSi2 06 was carried out on a 
large single crystal but we have no expectation for the 
ratio A\ : A2 because the details of the magnetic struc- 
ture are unclear. In Fig. U (b) we present the values of a 
derived from the asymmetry data for natural NaFeSi20g 
and the linewidth Ai associated with the oscillating com- 
ponent of the signal in the synthetic sample. The relax- 
ation rate a does not follow the same power law as the 
precession frequencies, suggesting that the static mag- 
netism does not emerge in as well defined a manner as in 
the synthetic samples, instead growing smoothly through 
the two transitions observed previously^ with no evi- 
dence for an intermediate collinearly ordered phase. We 
can estimate the range of magnetic fields at the muon 
sites using the relation Ai? = v^c/Tm ~ 0.6 T, which 
is larger than the value of ~ 0.2 T associated with the 
precession frequencies in the other two samples. To gain 
a further understanding of this we firstly calculated the 
dipole field distribution for two plausible model magnetic 
structures, ferromagnetic chains coupled antiferromag- 
netically, and antiferromagnetic chains coupled antiferro- 
magnetically. In both cases the moments were taken to 
lie along the chain direction. The results from these two 
calculations are very similar, with the muons sitting near 
the oxygen atoms linking the Fe octahedra and Si tetra- 
hedra, approximately a/4 from the Fe chains. Because of 
this we cannot distinguish the magnetic structure. The 
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FIG. 4: (Color online) (a) Muon oscillation frequencies, v 
(Eq.[2J, for LiFeSi2 0e and the synthetic sample of NaFeSi2C>6 
with fits to Eq.[3]described in the text. (Inset) The oscillation 
frequencies plotted against reduced temperature, t = (Tn — 
T)/Tn, showing the similarity of the trends approaching Tn. 
(b) The Gaussian relaxation rate, a (Eq. [4j| , for the natural 
NaFeSi2 0e and the linewidth Ai for the synthetic NaFeSi2 0e. 



second stage was to model the effect of local site dilution 
at the ~ 17 % of Fe 3+ sites which are not occupied by 
Fe 3+ ions. In this naive model we assume that all the 
dopants are non-magnetic and calculated the dipolc field 
of individual ions at muon sites. The distribution of fields 
is dominated by the effect of the closest iron moment to 
the muon site, which is around ~ 0.35 T, but averaging 
over the neighbouring sites leads to a distribution width 
AB ~ 0.25 T. Canting the moment direction towards the 
a-axis enhances this effect by up to a factor of two, but 
such a large canting is not consistent with the previous 
neutron diffraction results^ On this basis, site dilution 
alone cannot explain the magnetic field distribution ob- 
served in natural NaFeSi20g, but is consistent with the 
marginal change in the depolarization seen between the 
synthetic and natural samples. 



III. DISCUSSION 

Our heat capacity and /iSR results show that 
LiFeSi2 0g is a commensurate antiferromagnet in zero- 
field, in agreement with the neutron diffraction results 
reported previously^ The value of /? = 0.26(2) extracted 



from the temperature dependence of the oscillation fre- 
quencies suggests the magnetic ordering below Tn is in- 
termediate between two- and three-dimensional behavior. 
Synthetic NaFeSi2 0g shows similar heat capacity fea- 
tures to LiFeSi2 0g but the muon oscillations are heavily 
damped. This suggests a more complex magnetic order- 
ing where muons stopping at structurally equivalent sites 
experience a broad range of magnetic fields. The temper- 
ature dependence of the precession frequency in the syn- 
thetic NaFeSi206 sample is broadly consistent with that 
seen in the LiFeSi2 0g sample. This is seen more clearly 
when plotting the oscillation frequencies in the two sam- 
ples against reduced temperature, t — (Tn — T)/Tn, 
shown in the inset to Fig. H (a). Our heat capacity mea- 
surements on natural NaFeSi2 0g show that very little 
entropy is associated with the two previously identified 
magnetic transitions and instead short-ranged magnetic 
correlations build up in the quasi-one-dimensional chains 
over a temperature range extending well above Tn- The 
/iSR data appear similar to those for the synthetic sam- 
ple, except that the oscillations have become incoherent. 
This suggests that the impurities change the magnetic 
ground state, either by breaking up the intra-chain or- 
dering 3 ^ or by inducing a significant staggered magneti- 
zation around the impurity sites^i. Because of the com- 
plexity of this natural system it is not possible to sepa- 
rate these possibilities, nor be certain which of the two 
phases observed in the natural sample is present in the 
synthetic sample. From the shapes of the hump due to 
short-ranged order in each sample we can estimate^ in- 
trachain exchange constants ( J Ll ~ 7 K, J Na ~ 6.5 K) 
roughly consistent with the calculations of Streltsov and 
Khomskii 25 , though the comparison is complicated sig- 
nificantly by the interchain exchange. 



We can also compare our results to those reported on 
the other quasi-one-dimensional multiferroics LiCu202 12 
and Ca3(Co, Mn^Oe* 3 . The analogy with LiCu202 is 
somewhat closer, particularly for natural NaFeSi206, 
since there are two closely spaced magnetic transitions 
bounding a magnetically but not ferroelectrically or- 
dered intermediate phase. Our heat capacity measure- 
ments suggest that natural NaFeSi20g has consider- 
ably more one-dimensional magnetic interactions than 
LiCu202, but the magnetic structure is likely to be 
similar. Ca3(Co, Mn^Og has a similar hump in its 
heat capacity but below the onset of ferroelectricity 
and without pronounced features^ Both LiCu202 and 
Ca3(Co, Mn)2 06 show considerable evidence for disor- 
der influencing the multiferroic properties, due to Li 
non-stoichiometry and on-chain site disorder respectively. 
Comparison with the /uSR data for Ca3(Co, Mn^Oe 32 
shows that natural NaFeSi20g may indeed show some 
influence from incommensurate magnetism, on-chain site 
disorder, and correlations along the chains, as the muon 
precession we expect for the ordered states is evidently 
incoherent. 
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IV. CONCLUSIONS 

In conclusion, we have investigated the magnetic prop- 
erties of LiFeSi2 06 and both natural and synthetic 
NaFeSi2C>6 using muon spin relaxation and heat capacity 
measurements. LiFeSi20g enters a commensurate anti- 
ferromagnetic state below Tn = 18.5 K whereas both 
NaFeSi2C>6 samples appear to be incommensurate. In 
the natural sample, impurities within the Fe chains dis- 
turb this state and no coherent muon oscillations are ob- 
served. An unusual hysteresis is apparent in applied field 
heat capacity measurements of LiFeSi2 06 suggesting that 
magnetic correlations can be locked in well above Tn. 
Two magnetic transitions are apparent in the heat ca- 
pacity of the natural NaFeSi2C>6 but no related features 
occur in the /xSR data, which show a gradual build-up of 
static magnetism with decreasing temperature, or in the 
measurements on the synthetic sample, where only one 
transition is evident. Investigations of impurity effects 
in model one-dimensional chain magnets could shed con- 
siderable light on this enigmatic behavior. Further work 



would be worthwhile to make synthetic single crystals of 
NaFeSi2C>6 suitable for dielectric measurements to clarify 
whether the multiferroicity is indeed being triggered by 
disorder, and also to determine the magnetic structures 
of both synthetic and natural NaFeSi2C>6 to gain a bet- 
ter understanding of the multiferroic mechanism in this 
compound. 
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